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The proceeding comment raises a few points concerning our paper Dantchev et al, Phys. Rev. 

E. 89, 042116 (2014). In this reply we stress that while Refs. Diehl et al. EPL 100, 10004 (2012) 
and Phys. Rev. E. 89, 062123 (2014) use three different models to study the the Casimir force for 
the Ofn —>■ oo) model with free boundary conditions we study a single model over the entire range 
of temperatures, from above the bulk critical temperature, T c , to absolute temperatures down to 
T = 0. The use of a single model renders more transparent the crossover from effects dominated 
by critical fluctuations in the vicinity of the bulk transition temperature to effects controlled by 
Goldstone modes at low temperatures. Contrary to the assertion in the comment, we make no 
claim for the superiority of our model over any of those considered by Diehl et al. We also present 
additional evidence supporting our conclusion in Dantchev et al., Phys. Rev. E. 89, 042116 (2014) 
that the temperature range in which our low-temperature analytical expansion for the Casimir force 
increases as L grows and remains accurate for values of the ratio T/T c that become closer and closer 
to unity, while T remains well outside of the critical region. 


Before responding to the points raised in the com¬ 
ment by Diehl et al. [l] we summarize our article Q. 
Making use of the spherical model version of a lattice of 
fixed length spins with nearest neighbor interactions, we 
calculated the thermodynamic Casimir force of a three 
dimensional cubic lattice of such spins that is infinite 
in two dimensions and that has a finite extent, L, in 
the third. We applied free boundary conditions in the 
third dimension, which are more physically relevant than 
the more tractable, and well explored, option of peri¬ 
odic or anti periodic boundary conditions. Because of 
this choice it was necessary to implement separate spher¬ 
ical constraints layer by layer in that direction. Our 
calculations yield numerical results for the Casmir force 
that span all interesting temperature ranges, from above 
the bulk critical point down to absolute zero. When 
\T—T C \L is not too large—that is, in the finite size scaling 
regime—our results are consistent with finite size scaling, 
in that the Casimir force per unit area takes the form 
F Cas {T,L)/A = k B T c L~ 3 ’d((T - T C )L ), with D a univer¬ 
sal function of the scaling combination x = (T — T C )L. 
All our results in this region agree with those obtained 
by Diehl and co-workers [1, 0|- By contrast, outside of 
that regime, and most particularly at low temperatures, 
we discovered that the Casimir force is inconsistent with 
finite size scaling, in that there is a contribution that 
cannot be incorporated into the scaling form above. 

We supplemented our numerical results with an ex¬ 
pansion about the zero temperature limit of the Casimir 
force. Making use of this expansion we obtained an ana¬ 
lytical expression for the Casimir force that is asymptot¬ 
ically exact at T = 0, that is also exact for the leading 
order expansion about that limit. Furthermore, the ex¬ 
pression agrees with our numerical results to an excellent 
degree of accuracy over a wide range of temperatures, 
excluding of course the finite size scaling region. It is im¬ 
portant to note that the model we studied was viewed for 
some time as analytically intractable [1;, 45]. Nevertheless, 


we were able to extract analytical results at temperatures 
well below the critical one, and in addition, provide indi¬ 
cations that it may be possible to derive exact results at 
the critical point. All the steps in our investigation are 
laid out in the paper so that anyone interested in doing 
so can reproduce our results, extend or critically evaluate 
them, or make use of our methods in another context. 

By contrast, looking at the same geometry, Diehl et 
al. utilize the Ginzburg-Landau-Wilson (GLW) effec¬ 
tive Hamiltonian, studying two versions of the n —> oo 
Ofn) model, as described in the comment and in SI- 
They focus primarily on behavior in the finite size scaling 
regime, performing exhaustive and sophisticated analy¬ 
ses to extract results in that regime that are accurate to 
a remarkable number of significant figures. Taking limits 
in a particular way, they map on to earlier findings of 
Bray and Moore on the semi-infinite 0{n —> oo) model 
with Dirichlet boundary conditions. In the vicinity of the 
bulk critical point the authors focus on two GLW models, 
model A in which spins are confined to discrete layers per¬ 
pendicular to the bounding surfaces but are continuous 
distributed within each layer and model B in which spins 
lie on a two dimensional lattice within those layers. Un¬ 
surprisingly, corrections to scaling differ between the two 
models, while those models exhibit identical asymptotic 
critical behavior. For low temperature calculations they 
map onto a non-linear sigma model which for brevity we 
denote as model C. According to the authors model C, 
the analysis of which is not detailed, predicts asymptotic 
low temperature behavior that is consistent with finite 
size scaling, in that it can be expressed as a function of 
the scaling variable x. 

It is important to note that the three models above 
belong to the same universality class as the model we 
consider. It is also important to keep in mind that mem¬ 
bership in the same universality class does not guaran¬ 
tee identity of thermodynamic quantities outside of the 
critical regime. Even though a uniaxial antiferromag- 
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net and a gas-liquid system have the same equilibrium 
critical exponents at their respective critical points, a 
liquid freezes at sufficiently low temperatures and high 
pressures, while the typical antiferromagnet exhibits no 
corresponding behavior. No evidence has been presented 
that the model we consider has the same low tempera¬ 
ture thermodynamic Casimir force as the nonlinear sigma 
model, nor has any argument been presented that the 
nonlinear sigma model utilized in [3, hi will reproduce 
the low temperature Casimir force of either model A or 
model B. 

We stress that our study is performed on a single 
model. We make no claims for the superiority of this 
model over any of those considered in @, 0. Rather, 
we emphasize that our use of that one model to explore 
the Casimir force over the entire range of temperatures 
renders more transparent the crossover from effects dom¬ 
inated by critical fluctuations in the vicinity of the bulk 
transition temperature to effects controlled by Goldstone 
modes at low temperatures. 

As the comment lj makes clear, but which is not ev¬ 
ident in the articles [3, 0 , the low-temperature behavior 
studied there corresponds to T < T c but with —£ < 1, 
where their definition of £ is £ = —47r(i? — R c ) with 
R = f5J , which implies T/T c > 0.76, i.e., below but 
nevertheless still close to T c . In other words, the authors 
of the comment simply did not study “temperatures con¬ 
siderably below that of the bulk transition”, as asserted 
in0. 

We find that, if one fixes R ~ l/T large enough, then 
for any L the Casimir forces for different L’s cannot be 
made to fall on a single curve by invoking finite size scal¬ 
ing. This is in fact clearly shown in the plot displayed 
in (lj as well as in Figure 3 in our article 0. Further¬ 
more, we find that the numerically determined Casimir 
force is well reproduced by our analytical expressions. 
That is to say, those expressions accurately predict the 
Casimir for systems with fixed L at low temperatures. 
We agree with the statement in the comment that our 
“expansion never captures the correct asymptotic behav¬ 
ior”, i.e., the the scaling behavior of the Casimir force in 
the critical regime. On the other hand, the finite size 
scaling form obtained in 111 does not reproduce the 
correct true behavior of the Casimir force at fixed L for 
low enough temperatures. The absence of finite size scal¬ 
ing in the low temperature regime is by no means self- 
evident, given that when periodic or anti-periodic bound¬ 
ary conditions are imposed in the finite direction M 
finite size scaling holds down to T = 0. This leads to the 
infererence that the presence or absence of finite size scal¬ 
ing in the low-temperature, Goldstone mode dominated, 
regime strongly depends on boundary conditions. On the 
other hand, for all boundary conditions now studied the 
limit limi^oo lim T ^, 0 limA-i-oo[A 3 Fcas(r, A)/(fcsrA)] is 
universal; for periodic and Dirichlet boundary conditions 
this limit is given by the value of the scaling function of 
the force for the Gaussian model at its critical tempera¬ 
ture when subject to the those boundary conditions. 


By contrast, the results reported in i a , includ¬ 
ing those pertinent to their “low-temperature regime”, 
are formulated in terms of the scaling variable x = 
—4 tt(R — R C )L. We now understand that our original 
interpretation of that regime as encompassing tempera¬ 
tures close to absolute zero was in error. The error is un¬ 
derstandable given that it is natural to assume that a low 
temperature regime would refer to low absolute temper¬ 
atures, rather than temperatures for which T/T c > 0.76. 

With the above in mind we hope it is clear why we 
stated that we “partially corrected” the results of 1 in 
the low-temperature regime, while at the same time we 
are in full agreement with their results for the scaling 
regime. In fact, we make a special point in [2] of noting 
that our low temperature results are inconsistent with 
finite size scaling, so it is no criticism of those results 
to point out that they do not apply in the finite size 
scaling regime. In fact, we consider the inconsistency 
of our low temperature results with finite size scaling to 
constitute strong evidence that the Casimir force at low 
temperatures arises from different physical mechanisms 
(Goldstone modes) than does the critical Casimir force. 

With regard to the assertion in the comment that the 
results in 0] and 0 are “numerically exact” we note that 
0 contains a correction of the low temperature Casimir 
force as reported in (3); see the footnote on page 16 in Ref. 
0. As for the phrase “numerically exact,” the one inter¬ 
pretation that makes sense, as we see it, is that the au¬ 
thors utilize an algorithm that produces predictions that 
are as accurate as the implementation of that algorithm 
allows. By this interpretation our results are also numeri¬ 
cally exact. In any case, we did not dispute the numerical 
precision of the data presented in 1 and 0 for a fixed L. 
We simply pointed out that a related but not identical 
model analyzed at low temperatures provides informa¬ 
tion about the behavior of the thermodynamic Casimir 
force that is not subsumed in the analysis performed by 
the authors. Furthermore we are able to present exact 
closed-form analytical results for the Casimir force in the 
very low temperature regime (T asymptotically close to 
absolute zero). 

Our low temperature results are criticized in the com¬ 
ment, in which the argument is made that our reliance 
on an expansion in absolute temperature limits the ap¬ 
plicability of our expression to very low temperatures. 
However, it is a numerical fact that our analytical re¬ 
sults quite accurately reproduce the behavior of the 
Casimir force for reasonably large values of L and for 
all T < 0.8 T c , as shown in Fig. 1 of our article 0. In 
fact, the larger L is the better we find the approximation 
to be. We take this to be a strong empirical argument 
for the utility of the expansion. To further quantify our 
findings about the applicability of the low temperature 
expressions, the table below displays that values of T/T c 
below which the expression is accurate to both 5% and 
10%. The bounding ratios are shown for values of L 
ranging from 50 to 500. Note that when L = 500, 5% 
accuracy is achieved for all values of T/T c that are less 
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L 

50 

100 

150 

200 

300 

500 

5 % 

T/T c 

0.89 

0.901 

0.913 

0.922 

0.935 

0.949 

X* 

-19.6 

-34.8 

-45.5 

- 53.8 

-66.52 

-84.6 

10 % 

T/T c 

0.938 

0.947 

0.953 

0.958 

0.965 

0.974 

X* 

-10.5 

-17.8 

-23.4 

- 27.8 

-34.3 

-41.9 


than about 0.95. 

The accuracy of our expressions over a broad range of 
temperatures does lead one to ask why that is so. In Q 
we presented arguments that our analytical expressions 
should be accurate when — x = 4 tt(R — R C )L logL. 
As it turns out, the table above has x* increasing with 
L. That increase is slower than linear in L. The 5% 
data for x* are perfectly fitted by —21.98(1 + 3.06 logL — 
8.84 log log L). In any case, it is clear that as L increases 
our expressions remain accurate for values of the ratio 
T/T c that become closer and closer to unity. While the 
question of the exact range of applicability of our an¬ 


alytical expression is still open, the actual numbers do 
support our contention that the range of accuracy of our 
expression grows with increasing L. 

This leads us to the following statement in [x] : “our 
numerical results presented in Fig. 1 show that the range 
in temperature where the behavior is well described by 
the scaling function essentially does not depend on the 
thickness L, or equivalently, the range in the scaling vari¬ 
able x increases proportional to L". This is inconsistent 
with our claim for the expanding accuracy of our low 
temperature expression. Furthermore it conflicts with 
the standard implementation of finite-size scaling the¬ 
ory GM1, which is based on the expectation that the 
finite size scaling regime, in which the behavior of the 
finite system differs essentially from that of the infinite 
system, shrinks in absolute temperature as L increases. 
Our results, on the other hand, are in a full conformity 
with this expectation. 
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